Hyperfine structure of S states in Li and Be 
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A large-scale configuration-interaction (CI) calculation is reported for the hyperfine splitting of the 
2'^S and 3^S states of ^Li and ^Be^. The CI calculation based on the Dirac-Coulomb-Breit Hamil- 
tonian is supplemented with a separate treatment of the QED, nuclear-size, nuclear-magnetization 
distribution, and recoil corrections. The nonrelativistic limit of the CI results is in excellent agree- 
ment with variational calculations. The theoretical values obtained for the hyperfine splitting are 
complete to the relative order of and improve upon results of previous studies. 
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The hyperfine structure (hfs) of few-electron atoms 
has been an attractive subject of theoretical studies for 
decades, one of the reasons being a few ppm accuracy 
achieved in experiments on Li and Be"*" P, Inter- 
est in this topic was enhanced even further recently, in 
view of prospects of using hfs data to get an access to 
the neutron halo structure, the proton charge distribu- 
tion, and the nuclear vector polarizability, particularly 
for isotopes of Be"*" [1, 01 • Despite the considerable at- 
tention received, a high-precision theoretical description 
of hfs in few-electron atoms remains a difficult task. The 
main problem lies in the high singularity of the hfs inter- 
action and, as a consequence, in the dependence of the 
calculated results on the quality of the correlated wave 
function near the nucleus. 

Among numerous theoretical investigations performed 
previously for Li and Be+ , two apparently most accurate 
ones are the multiconfigurational Dirac-Fock (MCDF) 
calculation Q and the Hylleraas-type variational calcu- 
lation [3]. Both studies report good agreement with the 
experiment, but they are not entirely consistent with each 
other in treatment of individual corrections. The MCDF 
calculation does not include the binding QED effects and, 
in the case of Li, the nuclear magnetization distribution 
effect. The variational calculation yields accurate results 
for the nonrelativistic Fermi contact term but treats the 
relativistic effects in an effective way only, by rescaling 
the hydrogenic result. This indicates that neither of these 
studies is complete at the relative order of (a is the 
fine-structure constant). The aim of the present inves- 
tigation is to perform a high-precision calculation of the 
hfs splitting in Li and Be+, with a complete treatment 
of all corrections ~ . 

A possible way to accomplish this task would be to sup- 
plement the nonrelativistic calculation 6] with a rigorous 
evaluation of the relativistic correction, whose expression 
was recently derived by Pachucki Q . Such a calculation 
has not been performed so far. In the present work, the 
relativistic correction will be accounted for by means of 
the Dirac-Coulomb-Brcit Hamiltonian. 

The magnetic dipole hfs splitting of an energy level 
of an nS state is conveniently represented in terms of a 



dimensionless function Gn{Z) defined as 
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where /i is the nuclear magnetic moment, fiN — 
\e\/{2mp) is the nuclear magneton; m, m„, and M are 



nip, 

the masses of the electron, the proton, and the nucleus, 
respectively; / is the nuclear spin quantum number, and 
Z is the nuclear charge number. The function G defined 
in this way is unity for a non-relativistic point-nucleus 
H-like atom. 

Within the leading relativistic approximation, the elec- 
tron correlation can be described by the Dirac-Coulomb- 
Breit equation, which is solved by the configuration- 
interaction (CI) Dirac-Fock (DF) method in the present 
work. The many-electron wave function ^'(PJM) with 
the parity P, the momentum quantum number J, and 
the momentum projection M is represented as a sum of 
configuration-state functions (CSFs), 



*(PJi\/) = Cr^{-frPJM) . 



(2) 



The CSFs are obtained as linear combinations of the 
Slater determinants constructed from the positive-energy 
solutions of the Dirac equation with the frozen-core DF 
potential. The mixing coefficients Cr are determined by 
diagonalizing the Hamiltonian matrix. The hfs splitting 
is obtained as the expectation value of the hfs operator 
on the many-electron wave function The correspond- 
ing formulas are well-known, see, e.g., To perform 
a CI calculation, we devised a code, incorporating and 
adapting a number of existing packages 0] for setting 
up the CSFs, calculating angular-momentum coefficients, 
and diagonalizing the Hamiltonian matrix. The largest 
number of CSFs simultaneously handled was about a half 
of a million, with the number of nonzero elements in the 
Hamiltonian matrix of about 5 billions. A thorough opti- 
mization of the code was carried out, in order to keep the 
time and memory consumption of the calculation within 
reasonable limits. 

The dominant part of the hfs splitting in light atoms 
is delivered by the Dirac-Coulomb Hamiltonian. This 
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TABLE I: The Dirac-Coulomb-Breit part of the hfs splitting, 
in terms of G{Z). 





^max 


Li 2^5 


Be+ 


Coulomb 


1 


0.214 470 3 


0.390 159 9 




2 


0.2151678 


0.390 798 6 




3 


0.215 304 4 


0.390 938 7 




4 


0.215 346 2 


0.390 984 7 




5 


0.215 362 9 


0.391 003 8 




6 


0.215 3719 


0.391014 5 




7 


0.215 376 5 


0.391 020 2 




oo 


0.215 384 8(49) 


0.391030 4(61) 


Breit 




0.000 015 9 


0.000 038 6 


Total 




0.215 400 7(49) 


0.391069 0(61) 


MCDF [5] 




0.215 287 


0.390 984 


Hylleraas" [6] 




0.215 379(13) 


0.391 023(34) 



the sum of the nonrelativistic, the relativistic, and the 
nuclear-charge distribution terms. 



was the most demanding part of the calculation since a 
high relative precision was required. The one-electron 
orbitals for constructing CFSs were obtained by the 
dual-kinetic-balance (DKB) B-spline basis set method 
[lo| for the Dirac equation. For a given number of B- 
splines tIq, all eigenstates were taken with the energy 
< e < mc^(l + Zai?,nax) and the orbital quantum 
number I < ^max, where i^max was varied between 0.5 
and 6 and Zmax, between 1 and 7. Three main sets of one- 
electron orbitals were employed in the present work: (A) 
20s 20p 19d 19/ 185 ISh with = 44 and E^^^ = 3, (B) 
14s Up Ud 13/ 135 13h 12i 12k with Ua = 34 and £;,nax ^ 
0.5, and (C) 25s25p24d with = 54 and E^^ax = 6. 
Here, the notation, e.g., 20p means 20pi/2 20p3/2. Cal- 
culational results were first obtained with the set (A) 
and then corrected for contributions of the higher partial 
waves with the set (B) and for a more complete represen- 
tation of the Dirac spectrum with the set (C). The set of 
CSFs used in the calculation was obtained by taking all 
single, double, and triple excitations from the reference 
configuration with at least one electron orbital with I < 1 
present. The triple excitations that were left out in this 
way were found to yield a negligible contribution. In- 
clusion of the Breit interaction into the Dirac-Coulomb 
Hamiltonian yields only a small correction in the case 
of Li and Be"*". Since the effect is small, it is sufficient 
to use a much shorter basis set for its evaluation, which 
simplifies the computation greatly. 

The results of our CI calculation of the Dirac-Coulomb- 
Breit part of the ground-state hfs in ^Li and ^Be+ are 
presented in Table HI The Fermi model was employed 
for the nuclear-charge distribution, with the nuclear- 
charge radn ^ <r2>i/2= 2.431(28) fm for Li and 
<j.2^i/2= 2.518(11) fm for Be. The uncertainties spec- 
ified in the table include the estimated error due to the 
incompleteness of the basis and due to the finite nuclear 
size. The error of the Breit part was found to be negli- 
gible. Our results are in reasonable agreement with the 
nonrelativistic variational results but deviate signifi- 



cantly from the MCDF values [1] . The comparison leads 
us to a conclusion that the dominant part of the rel- 
ativistic correction can indeed be accounted for by an 
effective scaling of the hydrogenic results, as was argued 
in A complete evaluation of the relativistic correc- 
tion within the Za-cxpansion approach, however, has to 
be performed along the way paved in [7|, which has not 
been done yet. 

Comparison with the results of would become pos- 
sible on a much higher level of accuracy if we identified 
the nonrelativistic part of our CI results. Such an identi- 
fication was carried out by repeating our CI calculations 
for different values of a (namely, three values with ra- 
tios a' /a ~ 0.9, 1, and 1.1 were used). For each value of 
a, the finite nuclear-charge distribution (FNC) correction 
was evaluated separately and subtracted from the CI val- 
ues. The point-nucleus results thus obtained were fitted 
to a polynomial in a, assuming the absence of the linear 
term. In this way, the CI results with the physical value 
of a were separated into three parts: the nonrelativistic 
point-nucleus contribution, the relativistic point-nucleus 
correction, and the FNC correction. The numerical re- 
sults for them are listed in Table HTl 

The FNC correction was evaluated for both the hy- 
drogenic wave functions and the CI many-electron wave 
functions. In the latter case, a series of the CI calcula- 
tions with different values of the nuclear-charge radius R 
was performed and the FNC correction was extracted by 
a fit, using the analytical form of the R dependence [g]. 
It was found that, with an accuracy of ~ 0.5%, there was 
no screening effect on the relative value of this correction. 

The QED effects induce the largest correction to be 
added to the Dirac-Coulomb-Breit hfs value. For nS 
states of few-electron atoms, the QED correction can be 
written in the same form as for hydrogen [l3|. 



6Gn{Z) = :iGnZ)d + Za7T(\n2-^ 
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where G^^ is the nonrelativistic hfs value. The first 
three coefficients in the Za expansion ^ are the same 
as for hydrogen. The higher-order terms 021 and 020 
are different and not known at present. One can, how- 
ever, estimate them with their hydrogenic values [11, [13: 
a2i(2s) = -1.1675, a2o(2s) = 11.3522, a2i(3s) = 
-2.3754, and a2o(3s) = 9.7474. A 100% uncertainty 
is ascribed to this approximation. Essentially the same 
treatment of the QED correction was reported in ^ ; the 
QED results of 5] differ by ^^40% due to the neglect of 
the binding QED effects [i.e., the terms in ([3]) beyond the 
first one]. 

The nuclear structure effects have significant infiuence 
on hfs and should be taken into account. Their rigor- 
ous description is a demanding problem. The way for 
its solution was paved in recent studies [1, [l^ ■ Practi- 
cal realizations of this approach, however, are so far re- 
stricted by two- and three-nucleon systems and their 
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TABLE 11: Individual contributions to the hfs splitting, in terms of G(Z). The experimental values for the function G for Li 
were inferred from the original references by using the nuclear magnetic moment pi/^jv = 3.256 426 8(17) [l^ . 





■^Li 2^5" 


^Li S^S 


^Be+ 


9Be+ 3^5 


Noiir6lcit ivistic 


0.215 251 " 


168 340 " 


390 544 " 


335 066 " 


Ref. [6| 


0.215 254(4) 


0.168 351 (13) 


0.390 549 (9) 




Relativistic 


0.000 205" 


0.000 159° 


0.000 664" 


0.000 564" 


Finite nuclear charge 


-0.000 055'' 


-0.000 043° 


-0.000 139" 


-0.000119" 


Dirac-Coulomb-Breit 


0.215 401 (5) 


0.168 456 (9) 


0.391069(6) 


0.335 510(9) 


QED 


0.000 182 (4) 


0.000 143 (4) 


0.000 289(12) 


0.000 250(12) 


Bohr-Weisskopf 


-0.000 024(3) 


-0.000 019(2) 


-0.000 062 (17) 


-0.000 053(14) 


Specific mass shift 


0.000 002 


0.000 002 


0.000 002 


0.000 002 


Negative-continuum 


0.000 002 


0.000 002 


0.000 005 


0.000 005 


Total theory 


0.215 563(7) 


0.168 584(10) 


0.391 304 (22) 


0.335 714(21) 


Ref. [6] 


0.215 54(2) 


0.168 58 (2) 


0.39127(4) 




Experiment 


0.215 561 1 (1)'' 


0.168 60 (2)= 


0.391260(1)" 








0.171 5 (4)^* 


0.391 240 (6^^ 





" These three entries are inferred from the corresponding Dirac-Coulomb-Breit values- their sum is expected to be more 
accurate than each of the entries separately; Ref. [H; " Ref. [H; Ref. [H; Ref. 0] with ^l/^lN = -1.177 432(3) [l2|; 

=t Ref. with ii/iiN = -1.177 49(2) [11]. 



extension for more complex nuclei like ^Li and ^Be looks 
problematic. 

The most widely used approach up to now is to account 
for the nuclear magnetization distribution [the Bohr- 
Weisskopf (BW) effect] by means of the Zemach formula 
[l9| . which is simple and apparently model independent. 
Such approach ignores inelastic effects, which can yield 
a large contribution [l8| . and it is not clear what uncer- 
tainty should be ascribed to such results. In the present 
study, we calculate the BW correction within the single- 
particle (SP) nuclear model [IOjIMIj in which the nuclear 
magnetic moment is assumed to be induced by the odd 
nucleon. This model is expected to be reasonably ade- 
quate for ''Li since it reproduces well the observable nu- 
clear magnetic moment basing on just the free-nucleon g 
factors, the difference being only 15%. For ^Be, the devi- 
ation is four times larger and the SP approach is expected 
to yield worse results. 

Within the SP model, the BW effect can be accounted 
for by adding a multiplicative magnetization-distribution 
function to the standard point-dipole hfs interaction [2l[ . 
The distribution function is induced by the wave func- 
tion of the odd nucleon and is obtained by solving the 
Schrodingcr equation with the Woods-Saxon potential 
and an empirical spin-orbit interaction included. The 
parameters of the potential were taken from The 
BW correction was calculated for both the one-electron 
wave functions and the CI many-electron wave functions. 
It was found that, with a very good accuracy (< 0.5%), 
there was no screening effect on the relative value of this 
correction. Our calculational results are larger than the 
Zemach-formula values of i by ~ 10% in the case of Li 
and by ~ 30% in the case of Be. The Zemach-formula 
result of fs] for Be is larger than the one of Q by a factor 
of four, which is due, we believe, to a misinterpretation of 
the Zemach formula in Q . Our computational results for 
the BW correction are presented in Table [Til The error 



bars specified were obtained as the difference of the SP 
and the Zemach values and should be regarded as order- 
of-magnitude estimations of the error. We checked that 
similar evaluations of the nuclear effect on hfs in ^He"*" 
agree well with a much more elaborate calculation of [3] ■ 
The leading recoil contribution is given by the mass 
scaling factor {l-\-m/M)~'^ included into the definition of 
the function G in ([1]). The remaining correction (within 
the nonrelativistic approach) is due to the specific mass 
shift (SMS) and is very small for the S states. We calcu- 
late it by introducing the SMS term {m/M) J2i<jPi Pj 
into the Dirac-Coulomb Hamiltonian and taking the in- 
crement of the CI results with and without SMS {p is 
the momentum operator). Our results agree with the es- 
timates obtained in Q but are more accurate. For the 
2 and 3 ^5' states of Li, we obtain SG = 2.0(2) x 10^*^ 
and 1.9(2) x 10~^, respectively, which should be com- 
pared with 2(5) X 10-6 and 2(20) x lO'^ from [1], respec- 
tively. 

The negative-continuum contribution might be of some 
importance in calculations involving the operators that 
mix the upper and the lower components of the Dirac 
wave function. The hfs operator is of this kind, so we have 
to obtain an estimation for this correction. We calculate 
the negative-continuum contribution by employing the 
many-body perturbation theory to the first order. The 
same one-electron DKB basis set was used as in the CI 
calculations, with the only difference that all negative- 
energy eigenstates were taken. 

Our total theoretical values for the hfs splitting pre- 
sented in Table [ll] agree with the results by Yan et al. 
6} but are more accurate. The present theory agrees 
very well with experiments on Li, the only exception be- 
ing the experimental result which contradicts both 
the theory and the result of a more recent measurement 
. The comparison of theoretical calculations with the 
high-precision experiment for the ground state of Be"*" 
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is complicated by the existence of two different val- 
ues for the nuclear magnetic moment [T^ . The smaller 
value yields a better agreement with our theoretical re- 
sult, but there is still a 2a deviation present. Having in 
mind that the experimental results for the magnetic mo- 
ment are in significant disagreement with each other, one 
can surmise the presence of underestimated systematic 
effects in one or both of these measurements. We thus 
employ the comparison presented in Table |lT] to infer an 
independent value of the magnetic moment, which reads 
/i(^Bc)//iAr = —1.17730(6) and is somewhat smaller than 
the both values from [l^l . 

In summary, we have performed a large-scale CI cal- 
culation of the hfs splitting of the I^S and 3^5' states 
of Li and Be"*". The results obtained from the Dirac- 
Coulomb-Breit Hamiltonian agree with the previously 



reported nonrelativistic values but are more accurate 
due to a rigorous treatment of the relativistic correc- 
tion. The QED, nuclear magnetization distribution, re- 
coil, and negative-continuum corrections were evaluated 
separately and added to the Dirac-Coulomb-Breit value. 
Detailed comparison with the earlier calculations were 
made and some inconsistencies in their previous treat- 
ment of individual corrections were revealed. The cal- 
culational results for Li are in good agreement with the 
experimental data. For Be"*", the theoretical prediction 
deviates from the experimental value by 2 or 3(t, depend- 
ing on the value of the nuclear magnetic moment used. 
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